Voting and Catalytic Processes with Inhomogeneities 
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We consider the dynamics of tlie voter model and of the monomer-monomer catalytic process 
in the presence of many "competing" inhomogeneities and show, through exact calculations and 
numerical simulations, that their presence results in a nontrivial fluctuating steady state whose 
properties are studied and turn out to specifically depend on the dimensionality of the system, the 
strength of the inhomogeneities and their separating distances. In fact, in arbitrary dimensions, we 
obtain an exact (yet formal) expression of the order parameters (magnetization and concentration 
of adsorbed particles) in the presence of an arbitrary number n of inhomogeneities ("zealots" in 
the voter language) and formal similarities with suitable electrostatic systems are pointed out. In 
the nontrivial cases n = 1, 2, we explicitly compute the static and long-time properties of the order 
parameters and therefore capture the generic features of the systems. When n > 2, the problems are 
studied through numerical simulations. In one spatial dimension, we also compute the expressions 
of the stationary order parameters in the completely disordered case, where n is arbitrary large. 
Particular attention is paid to the spatial dependence of the stationary order parameters and formal 
connections with electrostatics. 
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I. INTRODUCTION 



Recently, much attention has been devoted to the field 
of nonequilibrium many-body stochastic processes 0. 
In particular the study of exact solutions of prototypi- 
cal models such as the voter model 0] has proved to be 
fruitful for understanding a broad class of nonequilibrium 
phenomena 0. In modeling nonequilibrium systems, it 
is often assumed that the underlying spatial structure 
is homogeneous. However, in real situations stochas- 
tic processes take place in the presence of imperfections 
(dislocations, defects, etc) that modify locally the inter- 
actions (see e.g. 0, and references therein). It is 
therefore highly desirable to take into account the effects 
of disorder, inhomogeneities and defects or other spatial 
constraints within simple and mathematically amenable 
models. Motivated by the above considerations, in a re- 
cent letter 0] , the properties of a paradigmatic nonequi- 
librium statistical mechanics model (the voter model) 
in the presence of one single inhomogeneity (a zealot) 
have been studied and it was shown that the presence of 
single zealot has dramatic effects on the dynamics and 
the steady state. For this model, in low dimensions, all 
of the agents eventually follow the zealot. Obviously, 
real systems are quite complex and the case of a single 
defect cannot be considered as being generic. To gain 
some insight on more realistic situations, we present here 
an approach allowing us to compute exact properties, 
in arbitrary dimensions, of a class of stochastic many- 
body systems in the presence of n competing inhomo- 
geneities. This study is carried out in the context of 
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two physically relevant systems: the voter model and the 
monomer-monomer catalytic reaction (in the reaction- 
controlled limit). We consider such a study as a fur- 
ther contribution toward the understanding of a class of 
disordered nonequilibrium many-body processes (where 
inhomogeneities would not be spatially fixed but would 
be randomly distributed). We show that the presence 
of "competing inhomogeneities" (in the sense of locally 
perturbing the otherwise homogeneous dynamics) gener- 
ally results in a space-dependent fiuctuating steady state. 
The amenable case where n = 2 is analytically studied in 
detail and the static and long-time properties of the order 
parameters are obtained and their spatial dependence are 
computed. The situation where n > 2 is investigated by 
numerical simulations. Also, in one spatial dimension, we 
are able to compute the stationary order parameters in 
the completely disordered case (i.e when n is arbitrary 
large). We therefore show how the stationary magne- 
tization/concentration depends on the dimensionality of 
the system, the strength of the inhomogeneities and their 
separating distances. In particular, we show that the lo- 
cal perturbation of the dynamics may give rise to subtle 
coarsening phenomena. In ID and 2D, when the density 
of the inhomogeneities is vanishing in the thermodynamic 
limit there is still coarsening in the system. Oppositely, 
when the density of the competing inhomogeneities is 
non-zero there is no coarsening, even in one and two di- 
mensions. We obtain an exact, yet formal, expression of 
the order parameters (magnetization and concentration 
of adsorbed particles) in arbitrary dimension. In dimen- 
sions d = 2, 3 we pay special attention to the radial and 
polar dependence of these quantities. Also, formal sim- 
ilarities with electrostatic systems are pointed out. The 
organization of this work is the following: In the next 
section we introduce the inhomogeneous voter model. In 
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Section III, we present the general mathematical set-up 
and the formal solution of the problem. In Section IV, 
we study analytically the voter model in the presence of 
two "competing zealots" in one, two and three dimen- 
sions and provide numerical results for the case where 
n > 2. In Section IV. B, for the one-dimensional case, we 
also derive the expression of the static magnetization in 
the completely disordered situation where n is arbitrary 
large. Section V is devoted to the study of the process 
of monomer-monomer catalysis reaction on an inhomoge- 
neous substrate, whose mathematical formulation is very 
close to that of the (inhomogeneous) voter model and in 
Section VI we summarize and present our conclusions. 



II. VOTER DYNAMICS IN THE PRESENCE OF 
COMPETING ZEALOTS 

The (homogeneous) voter model is an Ising-like model 
where a spin ("individual"), associated to a lattice site 
r, can have two different "opinions" <7r = ±1 0. The 
dynamics of such system is implemented by randomly 
choosing one spin and changing its state to the value of 
one of its randomly chosen nearest neighbors. In the 
(homogeneous) voter model, the global magnetization is 
conserved and the dynamics is Z2 symmetric (invariance 
under the global inversion cr^ —^ — Cr)- The importance 
of the voter model stems from the fact that it is one of a 
very few stochastic many-body systems that are solvable 
in any dimension. It is useful for describing the kinetics 
of catalytic reactions , for studying coarsening phe- 
nomena 0, 13 and also serves as a prototype model for 
opinion dynamics 0, 0, . 

Concepts inspired by statistical mechanics have al- 
ready been employed to some extent in the last two 
decades to mimic social issues Very recently vari- 

ants of the voter model and modern tools of noncquilib- 
rium statistical physics, such as various mean-field-likc 
approaches and exact methods 4,, .12, .15.] , numerical sim- 
ulations Hi 0, 0, 01 , formulation on random networks 
0,01 (see also references therein), were used intensively 
to quantitatively study further, both mathematically and 
numerically, collective phenomena, such as the opinion 
formation, inspired by socio-cultural situations. In this 
context, the voter model and its variants play a key role, 
as it is often used as a reference model. Despite of all 
these efforts, voter-like models have mainly be studied 
on homogeneous and/or translationally- invariant spatial 
structures. 

In contrast to most of the previous works, here we 
study, using exact analytical methods and numerical sim- 
ulations, a spatially inhomogeneous voter model. It is 
defined on a hypercubic lattice of size (2L -I- 1)'^, where 
individuals, labeled by a vector r having components 
—L<ri<L (with i — 1, . . . , rf), may interact accord- 
ing to the usual voter dynamics. In addition, we now 
consider that there are n zealots (labeled j = 1, . . . , n), 
occupying the sites {a^ = (aj, . . . , a^)}. These agents 



interact with their neighboring spins in a biased fashion. 
A zealot at site favors one of the opinions ej = ±1, 
i.e. it flips with an additional rate aj > (additional 
to the usual voter rate) toward his favorite state. As the 
zealots interact effectively with all of the spins on the 
lattice, there is a competition between them aiming at 
"convincing" as many spins as possible. Clearly, because 
the zealots perturb the dynamics locally, the system is 
disordered, not translationally invariant and the magne- 
tization is not conserved. 

According to the spin formulation of the model, the 
state of the system is described by the collection of all 
spins: S = {cr}- In this language, the dynamics of the 
model is governed by the usual voter model transition- 
rate 2, 5, 6] supplemented by local terms involving the 
zealots' reaction. The spin-flip rate, Wr{S) = •w{crr 
—(Tr), therefore reads: 



Wr{S) = 




Here the sum on right-hand side (r.h.s.) runs over the 2d 
nearest neighbors r' of site r and t = 1/d > defines the 
time scale. The probability distribution P{S, t) satisfies 
the master equation: 

^P(5,t) = J2 K(5'^)P(^^^) - Wr{S)P{S,t)] , (2) 

r 

where the state differs from configuration S only by 
the spin-flip of cr^. Using the master equation |(2Jl, in the 
thermodynamic limit i — s- cx), the equation of motion 
of the local magnetization at site r, denoted by Sr{t) = 
J2s '^r P{S, t), reads: 

-Srit) = Ar.Sr{t)+J2^^ fe- - ^a.(i)) '^..a. • (3) 

Here A,, denotes the discrete Laplace operator: 
ArSr{t) = -2dSr{t) +Y.r'^r'{t). Wc cau immedi- 
ately notice from (j^J that the stationary magnetiza- 
tion obeys a discrete Poisson-like equation: Aj.Sr{oo) = 
J2]=i'^j ('5'a3(oo) — ej) 6r a,i- There is an obvious and 
striking resemblance between this equation and the well- 
known equation for the electrostatic potential generated 
by n classical point charges located at {a-'}. Therefore, 
one may be tempted to formally identify Sr{oo) with an 
electrostatic potential and think that the problem could 
be solved easily. In fact, the problem is much harder since 
the quantities playing the role of charges depend them- 
selves on the magnetization. In other words, the problem 
of finding the stationary magnetization is isomorphic to 
the problem of determining the electrostatic potential in 
a discrete system where the value of the charges depends 
on the potential itself. Because of this fact, the calcula- 
tion of Sr {oo) cannot be inferred directly from the results 
known from electrostatics and the computations have to 
be carried out in a self-consistent manner, as described 
hereafter. 
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III. GENERAL SET-UP AND FORMAL 
SOLUTION 

In this section, we show how to compute the magneti- 
zation of the voter model in the presence of an arbitrary 
number of inhomogeneities (competing zealots) and pro- 
vide a "formal" solution of Eq. © • 

For further use, we introduce the following quantity: 
lis) = J^dt e-^*[e-2'^%^(2t).../,,(2t)] = L^is), 
where /„(2t) = /_„(2t) = ^ cos (qn) e^*'^"'^'? is the 
usual modified Bessel function of first kind 20]. The 
quantity Ir{s) can be rewritten in terms of Watson inte- 
grals, or "lattice Green-functions": 



Iris) = I-r{s) = 



. (27r)'^5 + 2[d-EtiCOsg,] 



, (4) 



where q = (qi,...,qd) is a d— dimensional vector. We 
also introduce the Fourier transform of the magnetization 



(5) 



Fourier transforming (|3Jl, we obtain the following equa- 
tion: 

j^Sqit) - -2d (^l-^Y.co8q}j Sgit) 



i=i 



(6) 



where A' (t) = aj (e^ — S'qj (i)) . Laplace-transforming 
Eq. (jnj, we obtain the following expression for the 
Laplace-Fourier transform of the magnetization: 



s + 2d|l - -5 E^=i cos gij 



(7) 



where A^{s) = dte~^*A>{t). For technical simplicity, 
we have considered that the system is initially in a state 
with zero magnetization: 5*^(0) = 0. Inverse Fourier 
transforming Eq. Q, we get the Laplace transform Sr{s) 
of the magnetization: 



Sr{s) = 



A'^{s) e 



i{a —r).q 



+ 2c!{l-iEtiCosg.} 



(8) 



As both right and left hand-side (l.h.s.) still depend 
on the Laplace transform of the magnetization (through 
A^{s) on the l.h.s.), to obtain an explicit expression for 
Sao (s)i we have to find a self-consistent solution of Eq.® 
for all of the a^ 's by plugging r = into Eq.ljSJ. Solv- 
ing the resulting linear system, in thermodynamic limit 
{L oo) we obtain: 



Sa.{s) 



E 



d'^q 



A'^{s)e 



'.{a —a-').q 



^ (2'^)^ + 2d{l-iEtiCOS(?.} 



, (9) 



which can be rewritten Yli i^^j-i + 17^) A^{s) — ej/s, 
where the symmetric nxn matrix Ad is defined as follows: 



Mji{s) 



TT id 



d'^q 



^ (2^)S + 2d{l-iEtiCosg.} 

4.-a'(s)=/a^-a.(s). (10) 



To obtain the two last equalities, we used the integral 
representation Q . We now introduce another symmetric 
n X n matrix, Af, defined by: 



M,A-^,{a})^M,As)+^-^, 



(11) 



and from it, using Eq.Q, one obtains A^ and Sai- 

= iE^^[-^"'(«,W)k. (12) 

e 

Sa.is) - I (^6,-^^Q[AA-i(.,{a})],,,^ (13) 

At this point, we can get an explicit expression for the 
Laplace transform of the magnetization by plugging back 
(|12|l into ||HJ). In the thermodynamic limit [L ^ oo), we 
have: 

^r{s) - ^ E ^iL.-ris)W-\s, {a})],M (14) 

and therefore, formally the magnetization is obtained by 
Laplace- inverting Eq. (|14(l : 



Srit) = 



1 

27ri 



X / ^^e^*5]Q4._.(s)[AA-i(s,M)],,,.(15) 

J C — iOD ^ j g 

This expression means that we have recast the problem 
of solving the inhomogeneous voter model in the presence 
of arbitrary many inhomogeneities into a well-defined lin- 
ear algebra problem whose main, but nontrivial, analytic 
difficulty resides in the inversion of the matrix J\f. The 
steady state of the magnetization for L oo can be 
directly inferred from Ea. (|14|l as follows: 

S'r(oo) = lims^o eeia:>-r{s)[Af^^{s,{a})]j^i. (16) 

The exact expression for the long-time magnetization is 
obtained by Laplace-inverting the s — > expansion of 
Ea. ((TH) . after having subtracted the static contribution 
Sr{oo)/s, and by paying due attention to the situations 
where the integrals are divergent. It is also worth 
mentioning that the properties of the modified Bessel 
functions of the first kind, Ir{t) 20], allow us to write a 
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formal and implicit solution of Eq. ^ ioi L 
reads : 



oo, which 



Sr{t) = ^^fc(0)n [e-2*4^_^^(2t)] 

k 1=1 

+ J2 /* dt'A^it - t') n i2t')] .(17) 

j i=l 

To solve it explicitly for Sr{t), one has to Laplace trans- 
form (|17|l and then solve the resulting linear system 0|, 
which is equivalent to the procedure described above. 
The expression (|17|l is advantageous if one is interested 
in the global magnetization of the system. In fact, as we 
consider an initially homogeneous and "neutral" system 
{Sk{0) = 0), using Eqs H17|) . the global magnetization of 
the system can be written: 

M{t) ^ E^-^w^E / ^^^'M 

k 3 = 1 

= ^a, / dT [e,.-^„,(r)], (18) 

where we use the identity YlV=-oa ^kif) = e* 0|. 

The situation considered here is particularly interest- 
ing when the zealots favor different opinions and there 
is an effective competition occurring in the system. In 



this case we expect nontrivial nonequilibrium space- 
dependent steady states. Of course, we can easily check 
that in the presence of one single zealot [n — 1) located 
at site 0, with strength ai = a and ei = 1, we recover the 
results reported in Reference In this case we simply 



have: 



a/o(s) 



1 



we recover Sr{s) = 



and, together with ||SJ), 
In Ref . 4] , one of us has 



s(aio[s) + l) ■ 

shown that in low dimensions the voter model with only 
one zealot evolves toward the unanimous state favored 
by the inhomogeneity. 



IV. THE VOTER MODEL IN THE PRESENCE 
OF TWO COMPETING ZEALOTS 



In this section we specifically consider the case where 
two competing zealots are present (j = 1,2): One, with 
strength ai = a, located at site = Q and the other 
located at site a? ^ x with a strength a2 = /?. This 
case is explicitly tractable and displays interesting fea- 
tures, which turns out to be generic for the case n > 1 
as illustrated by numerical simulations. For this case, 
7o(s)-ha-i 4(s) 

Lis) io{s) + i3- 

and therefore, using Ea. (|14|l . we infer the expression of 
the Laplace transform of the magnetization at site r : 



we have J\f = 



for L 



oo. 



J 



Sris) 



r(s) + P<^2lr-xis) + |-^r(s) (el^(s) " e24(s)) + ^r- x(s) (e2^(s) - ei4(s)) I 
s[l + {a + /3)/o(s) -f «/3(/2(s) - /2(s))] ' 



(19) 



where ei^2 = ±1- Obviously, the inhomogeneous system 
with two zealots is interesting in the case when ei = — e2. 
In fact, it is clear from Ref. Q that in ID and 2D the 
condition ei = £2 implies that Sr{oo) — ei- In this sit- 
uation, the long-time approach toward the unanimous 
steady state is Srit) — Sr{oo) ~ At^^^^ in one dimension 
and Sr(t) — Sr{oo) ~ B/hit in two dimensions. Thus, 
in low dimensions, when ei = £2, only the long-time am- 



plitudes A and B change with respect to the case where 
n = 1 and £ = £1 P|. 

From now on, without loss of generality, we consider 
the more interesting situation when there is a competi- 
tion between the zealots: £1 — —£2 = 1. Namely, the 
zealot at the origin favors the +1 opinion, whereas the 
zealot at site x favors the opposite —1 state. In this case, 
Eq. (|19|l simplifies as follows: 



g ^ air is) - Pir-.js) + aP iUs) - ir-^{s))iio{s) + Ljs)) 

s[l + (a + /3)/o(s)+a/3(/2(s)-/2(s))] 

I 

Different questions can be asked here: What is the zealots ? How does the opinion of a randomly picked spin 
range of influence of each zealot ? How "efficient" are the 
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FIG. 1: (a) Graphical representation of a microscopic config- 
uration of the spins on a one-dimensional chain. The zealot 
favoring the -1-1 opinion with a strength a, indicated by a dot 
and a larger up-spin, is at the origin. On the right of the 
origin, at a distance x, the other zealot, indicated by a square 
and a larger down-spin, favors the —1 state with a strength /3. 
(b) Typical ID stationary magnetization profile 5*^(00) (de- 
noted simply Sr in the figure) versus r in the thermodynamic 
limit. On the left of the origin and the right of the other 
zealot, the static magnetization reaches two plateaus with 
heights given by Eqs. 1231 and (1241 . Between the zealots, 
the stationary magnetization varies linearly according to Eq. 

evolve with the time, and what will be its final opinion ? 
These questions will be answered in the next sections by 



explicit calculation of the stationary magnetization and 
its long-time behavior. 



A. Results in ID 

First we focus on the one-dimensional situation and 
consider the case when both competing zealots are sep- 
arated by a finite distance x [See Fig.nja)]. It is worth 
studying the properties of the one-dimensional version of 
the inhomogeneous voter model because of its physical 
implication for the catalysis (see Section V) and its close 
relationship with the Ising model with Glauber dynamics, 
which is an important theoretical model, known to have 
many physical applications 0,3. In fact, in the absence 
of zealots the one-dimensional voter model coincides with 
the Glauber-Ising model with zero temperature dynamics 

In ID, one computes explicitly /^(s) in Eq. as 
follows ,2a,,21J: 

Iris) = Iris) = ^ , , , -, (21) 

where r — \r\. We see that in ID Iris) diverges for small 
s as s^-'-/^. 

Without loss of generality we consider the situation 
illustrated in Fig. and thus, from Eqs. (j^OJ), the 
long-time expression for Srit) in the case where r e [0, x] 
is 



Srit) 



a - /? - a/3(2r - x) 
a + P + al3x 



(a + /3 + af3x)\/TTt 



2{a -13- aPi2r - x)} + afix {/3(a; - r) - ra] 



a + /? + a(3x 



For the spins on the right of the origin, with a; < r < 00, we find 



Srit) 



1 



a — (3 — a/Sx 
a + j3 + aj3x) (a + /3 + a/3a;)v^ 



2(0; — /3 — a(3x) — a^(3x^ 
a + P + a(3x 



+ ar -\- Pix — r) 



whereas for the spins on the left of the origin, with < r < 00, we find: 



S-rit) 



a — (3 + a[3x 
a + l3 + aj3xj (a + /3 + a^x) ^Trf 



2(a - /3 -t- al3x) + afS'^x'^ 



a + P + a(3x 



+ ar + /3(r — x) 



(23) 



+ ar — /3(r + x) 



(22) 



(24) 



Finally, when both r ^ 00 and t — > 00, Iris) — > 
e^"^^^ i2^/s) . Using this expression in Ea. (|2U|l . as in 
Ref.[3, we obtain a scaling expression for the magne- 
tization: 



S±rit) 



a — [3^ a/Sx 
a + l3 + a/Sx 



erfc 



(2^/^) 



(25) 



where erfc(a;) = 2 



dy 



is the usual complemen- 



tary error function. We infer from H22(l that in the 
finite interval separating the two zealots, the station- 
ary magnetization profile decays linearly with a slope 
— 2a/3/(a + f3 + a(3x). Outside from this interval, the 
final magnetization is uniform on the right and left 
hand side from both inhomogeneities. In fact, ^2'A\ and 
(|24|l show that the static magnetization of the spins is 
^±-(00) = frffiff (see Figs.[ll;b),|2l. These plateaus 
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FIG. 2: The stationary distribution S'r(oo) on a L = 1024 
lattice with two competing zealots. The zealot favoring the 
positive opinion (dot) is located at the origin with a — 0.02 
and the other one favoring the negative opinion (square) is at 
r = 430 with (3 = 0.01. The agreement with the theoretical 
results for an infinite system is excellent. 



differ significantly from the values =Fl only when when 
the product a/? is comparable to x~^. Therefore in ID, 
the final stationary solution, which is summarized on Fig. 
^b), is polarized and can be understood as being the 
solution of a discrete one-dimensional electrostatic Pois- 
son equation with peculiar boundary conditions. In fact, 
it is well known that in ID the electrostatic potential 
varies linearly with the distance to the charges. Here, 
the nontrivial part of the analysis is to compute in a 
self-consistent manner the heights of the plateaus. All 
these profiles are approached algebraically in time, i.e. 
Sr{t) — Sr{oo) — At'^l"^ (as in the case with only one 
zealot 0), with amplitude depending nontrivially of all 
parameters of the system A = A{a,P,x)r. Obviously, 
because there is a distance x separating the zealot at the 
origin from the other, the expression for Srit) is not sym- 
metric with respect to the site 0. We can notice that the 
expressions (I22|I . (I23|) and (|24|l simplify when the strength 
of the zealot is infinite {a = (3 = oo). In this case, the 
zealots have a final magnetization 5o(oo) = —5*2,(00) = 1. 

Result tells us that for spins infinitely far away 
from the zealots, the magnetization evolves as a smooth 
scaling function of the variable u = This scaling 

function differs from zero (the initial condition) after a 
long time {i.e. t ^ ^ 00), when the variable u has 
a finite value. It follows from Eqs.|j2Sl),lElI) and 
that in ID the effect of the zealots is felt and propa- 
gates as i^/'^ 00. For large time and distance, when 
1 <C t ^ r^, we see from Ea. (|^ that Sr{t) is still close 
to its initial value. When t ^ r"^, all the agents approach 
as the active and fluctuating stationary magnetiza- 
tion (123 • From Eqs. and ((TH|l . we can infer the 
long-time behavior of the global magnetization in the 
system. As a(l - 5o(t)) - /3(1 + 5,(i)) c 



when a ^ (3, the average number of voters following the 
strongest zealot evolves (at long-time) as the square-root 

of time: M{t] ~ l^aT^a \ -■ This result implies that 
the time T necessary for the strongest zealot to domi- 
nate (on average) the whole ID system scales as T ~ L^, 
where L — > 00. When a — (3^ the system is exactly sym- 
metric with respect to x/2, and in average there are as 
many -f 1 spins than —1 ones in the whole system. 

On Fig.|21we show the stationary magnetization Sr{oo) 
on a finite lattice with L — 1024 for two competing 
zealots obtained from Monte Carlo simulations. For sim- 
ulating the model we use random sequential dynamics 
by picking randomly an "active" site ( either one of the 
zealots or a site that has at least one nearest neighbor 
in a different state) and flipping it with a rate given by 
Eq. J^l- The time after an attempt for a flip is updated 
with the amount 1/Na, where Na is the number of ac- 
tive sites before the current update. To account for the 
fact that the simulations are on a finite lattice, where the 
spin at the left (right) boundary site has only one near- 
est neighbor on the right (left), the spin-flip rate at the 
boundaries is modifled such that it depends only on the 
state of a single neighbor. The first 2 x 10® Monte Carlo 
steps (MCS) are discarded and typically we sample the 
configurations on the lattice every 5000 MCS for the next 
5 X 10^ MCS. The stationary distribution for S'r(oo) ob- 
tained from the simulations is in an excellent agreement 
with the theoretical values obtained for a inflnite lattice 
and sketched on Fig.QIb). 

Fig. 131 shows the result from Monte Carlo simulations 
on a relatively small (L — 8192) lattice for various av- 
erage quantities. The long time behavior of the local 
magnetization (5S'o(i) = 150(00) — So{t) and SSx{t) = 
Sx{oo) — Sx{t) clearly show the long time behavior, 

in agreement with Ea. H22() . In Fig. |3]we also report nu- 
merical results for the average number of interfaces {i.e. 
two neighboring sites with different spins) vs. time. This 
quantity gives us a good qualitative and quantitative pic- 
ture of the coarsening of the system. Fig. El shows that 
the average value of the interfaces, which equals to the 
number of the clusters of -1-1 and —1 spins, evolves as 
t~^/^ before saturating at a small non-zero value. One 
can notice that for a long time the system evolves and 
coarsens as in the homogeneous voter model Q , but due 
to the presence of the two competing zealots, subtleties 
appear at long times. In fact, one has to distinguish be- 
tween the three possible situations for the coarsening: (i) 
when we have n < 2 (i.e. none or only one zealot on the 
lattice) , there is the usual coarsening (an infinite domain 
spans the entire system) Q ; (ii) when 2 < n and the den- 
sity {n/L) of the competing inhomogeneities is zero for 
L 00, there is still coarsening in the sense that the size 
of the different domains formed increases with the size of 
the lattice but never spans the entire lattice; (iii) when 
the density of the competing zealots has a non-zero value 
in the thermodynamic limit, there is no longer coarsen- 
ing as the formation of large domains is prevented by 
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FIG. 3: Coarsening on the one-dimensional model with two 
competing zealots. The figure shows the average number of in- 
terfaces vs. time, the average magnetization of the two zealots 
(see the text) So{t) and Sx{t), and also 5So{t) = So(oo)— S'o(t) 
and SSx{t) = Sx{oo) — Sx{t)- The simulation is on L = 8192 
lattice for a = 0.5, /3 = 0.2 and x = 3000 and the continu- 
ous lines shown have a slope —0.5, as predicted by Eq. 11221 . 
For this choice of the parameters, the average number of in- 
terfaces decays algebraically toward a small but finite value 
(here, « 2.0 x lO"*). 



the interaction with the numerous (competing) inhomo- 
geneities. 

After having discussed in detail the case n = 2, we 
would hke to point out that in one spatial dimension it 
is possible to compute the stationary magnetization for 
an arbitrary number n of zealots in a more direct and 
intuitive fashion than relying on Ea. H16() . In fact, let 
us consider that the zealots, labeled by j — 1, . . . , n are 
at sites — oo < ai < a2 < • • • < a„ < oo. By plugging 
the ansatz that the stationary magnetization between the 
sites and a^"*"-^ reads Sr {oo) = SaJ (oo) + ^j{r — a^) 
into ArS'r(oo) = — X]J=i o^i (^j ~ "^a^ (oo)) (5^ , where 



we have introduced 7^ 



S^j + i (oo)-S^j (00) 



and Xi 



we obtain: 



ll5r,ai + (72 - ll)6r,a2 H 1" {ln-1 - 7n-2 )'5r,a„_ i 

n 

ln5r,a„ = ^ OLj {SaJ (oo) - Ej) S^^ai ■ (26) 




FIG. 4: An example for numerical simulation of the case with 
4 zealots on a L = 1024 lattice (see the text). The bias of the 
zealots from left to right is 0.01 (negative), 0.02 (positive), 
0.013 (negative) and 0.003 (positive). 



netization at each sites < a' < a": 

Sai{oo) 



£2 



ai 

72 - 71 



S'a"-i(oo) 
5'a"(oo) 



a2 



In-l — 7n-2 



7n-l 



(27) 



Of course, in each of these equations for Sai {oo), 
the right-hand-side depends on Sai {oo) and 5'aj+i(oo) 
through 7j . The equations H27I) are therefore a set of cou- 
pled linear equations that can be rewritten as PS — v, 
where P is a nxn band matrix, which only non- vanishing 
entries are 



p 


= -i^J 




I < j <n 




= Xj, 


1 < j < 71 






~ ^i~i 


I < j <n 




Phi 


= -(1- 


\- aixi) 




P 

^ n.n 


- -(1- 


OLn-^n— 1 ) 




Pl,2 


— P 

— n,n- 


-1 = 1, 


{ 



(28) 

and S and v are column vectors which components are 
respectively 

S.J = Sai{co), l<j<n 
vi — —tiaiXi 

-e,a,x,_i a;, , 1 < j <n 



-CnO- jiXji . 



(29) 



Solving these equations, we obtain the stationary mag- 



Therefore, the solution of (|27|) is obtained by inverting 
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the band matrix P and reads: 

n 



(30) 



k=l 



Having solved (at least formally) the set of equations H27|) 
giving the stationary magnetization at each site a-', the 
general one-dimensional stationary magnetization in the 
presence of n zealots simply reads: 

• If r < a^: 

Sr{oo) = Saiioo). (31) 

• lia^ <r < a^+i (1 < j < n): 

Of \ Q f \ : -Sai + i (oo) - S'ai (OO) - 

5*^(00) = S'^^oo) + _ ^ (r-a^). (32) 

• If r > a": 

S',.(oo) = S'a"(oo). (33) 

As an example, let us consider the case where there are 
four zealots on the chain, as illustrated in Fig. 0| This 
figure shows that the one-dimensional stationary magne- 
tization profile is a piecewise function, as predicted by 
Eqs. 1211-1221. When n = 4, as in Fig. H Eqs. (EJ ex- 
plicitly read : 

5*02(00) - (1 + aia;i)S'ai(oo) = -eiaia;i 

XiSa^ioo) - {Xi+ X2+ a2XiX2)Sa2{cC,) + X2Sa^{co) 

= —€20^2X1X2 

X2Sai{<X>) - {X2+ Xz+ azX2Xi)Sai{<X>) + X:iSa-i{<x) 

= -6303X2X3 

Sa3{oo) - {I + a4X3)Sai{oo) = -6404X3 (34) 

The set of Eas. (j34(l can be solved explicitly and gives 
rise to very cumbersome expressions. Plugging into the 
latter the values corresponding to the system simulated 
in Fig. H i.e. ai = 0.01, ei = -1, a2 = 0.02,62 = +1, 
tta = 0.013,63 = —1 and a4 — 0.003,64 = +1, and 
xi — 230, a;2 — 240, = 130, we obtain: 5ai(oo) — 
-0.529,5^2(00) = -hO.556, 5^3(00) = -0.441, S-^i (00) = 
—0.0367. These values can be compared to the re- 
sults of the simulations, reported in Fig. ^ , where we 
obtained S'ai(oo) = -0.53 ± 0.01, S'a2 (00) = +0.55 ± 
0.01, 5^3(00) = -0.45 ± 0.01,5^1(00) = -0.04 ± 0.005. 
These comparisons show that there is an excellent agree- 
ment between the theoretical values predicted by the so- 
lution H30|l of the system H34() and the numerical results. 
This agreement is somewhat surprising as the simula- 
tions reported in Fig.^have been carried on a relatively 
small system (L — 1024), whereas all the theoretical re- 
sults (|27|) - (|34|l have been derived in the thermodynamic 
limit. This fact indicates that our analytic results may be 
quantitatively accurate even for large, but non-infinite, 
systems. In the limit where the strength of the zealots 



is Q!i = • • • = an ~ 00, all the expressions simplify 
and it folows from 127|1 that Sai {00) = ej, while, for 

<r< a^+\ 5,(c») = 6, + (^^) (r - a,). When 

ai = • • • = a„ = 00, this ID system can be related to 
the one-dimensional spin model with Glauber dynamics 
(at zero-temperature) in the presence of quenched ran- 
dom fields of infinite strength pla |: in the voter language, 
the situation considered by the authors of Ref . would 
correspond to the case where at each site j a "voter" 
would have a probability p to be a zealot favoring the 
opinion ej = ±1 with strength aj — 00 and would have a 
probability 1 — 2p to be an ordinary agent. The (slight) 
difference between such a model and the one studied in 
Ref. 10 is the fact that each zealot (even when he is en- 
dowed with an infinite strength) can be "forced" to flip 
by his two neighbors, while in Ref.T^ the (random) mag- 
netic fields pin the spins along their direction. However, 
as aj — 00, each zealot j rapidly flips back to his prefer- 
able opinion ej and thus both models are very close and 
display the same stationary magnetization. 

We also would like to emphasize that the results H27|l . 
(|30|l provide the exact magnetization of the completely 
disordered one-dimensional voter-model, where each site 
is endowed with a speciflc spin-flip rate. In this case, 
one would have n = L 00 zealots in the system with 
Xj = a^+^ — = 1, and the structure of the matrix P is 
rather simple [see Ea. H28l) ]. 



B. Results in 2D 

In two dimensions, the integral of Eq.l^J is also diver- 
gent in the long-time regime s — > and therefore its main 
contribution arises from = qf + q2 0. In this sense, 
we first expand Eq.l^J for small s in the case when r = 0: 



» -— Ins, 



(35) 



More generally, for r :s> 1, we have (see Ref.0) 
Ir{s) > where Ko{x) is the usual mod- 

s— >0 

ified Bessel function of the third kind |20|. Using the 
small argument expansion of such a Bessel function we 
find that the long-time behavior for t ^ ^ 1 is given 
by 



iris) 



1 



r^^O 47r 



[ln(r2s) + 2{7-ln2}] , (36) 



where 7 = 0.5772156649 . . . denotes the usual Euler- 
Mascheroni's constant. From the expression (|20|l . when x 
is sufficiently large to use Eq. H36|) , we obtain the station- 
ary magnetization of the zealots: 6*0(00) ~ " j^j 

and Sx{oo) — Jfi |" ^ . Interestingly these expres- 



sions resemble to the ones obtained in ID [see Eqs H23|l . 
(|24|l ] . The only change is in the dependence on sepa- 
rating distance: With respect to the ID case, one has 
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Sr>0 



Sr>0 




FIG. 5: Sketch of the typical 2D spatial dependence of the 
stationary magnetization when L oo. At the origin, in- 
dicated by a dot, is the zealot favoring the state +1 with a 
strength a = 1. At a distance a; S> 1, indicated by a square, 
is the zealot favoring the state —1 with a strength /3 ~ 0.9. 
According to Eq. 113711 . the agents within the disk of center 
c ~ 2a; and of radius R ~ 1.42; have a negative final magne- 
tization (denoted simply 5*^ in the figure). Outside the disk, 
the final magnetization of the agents is positive (see the text), 
while on the circle the agents are in a "neutral" final state. 
The static magnetization ^^(oo) exhibits both radial and po- 
lar dependence. 



X ^ ^^nx. When r 1 and |r — a;| ^ 1, from H20() . us- 
ing Eas. l|35|) and the stationary magnetization reads 
(see Fig. 5): 



Sr{oo) 



/J _ M In 

' TT |r — cc| 



r>i,|r-x|>i a /? + ^ [Inx + 7r(7 - ln2)] 



t(37) 



Far away from both zealots, and in the case of sufR- 

" ^ X ^ 1, this expression 



far away irom ooin zeaio 
ciently separated zealots, i.e. 
simplifies: 



Q/3 : 



cosf 



a + (3+^[\nx + n{-r-\n2)y 



(38) 



where cos 61 = — . 

rx 
X cos 



ln(r/|r — x\) 



Here we used the fact that 
-I- 0{{x/r)'^), when r > a: > 1. 



These results show that, because of the competition 
between the two zealots, the stationary magnetization 
is a fluctuating steady state exhibiting nontrivial ra- 
dial and polar dependence. Also, when a = (3 — oo, 

__ In 

Eq.(E2l) reduces to 5,(oo) — — > inx+^(;-in2) and 

5o(oo) = -Sxioo) = 1. 

Regarding the dynamical behavior, in the regime 
where t ^ max(|r — a;p,r^), the long-time behavior of 
the magnetization is the following: 



Sr{t) — 5*^(00) 



1 

hii 



In 



^ lnp^{ln(x/2)+7} + 2(a-/3)(7-ln2)' 



Q/3 



{Inx + 7r(7 - In 2)} 



(39) 



> 1, the 
approach to- 

— Sr{oo) ~ 

For t > r^. 



In the situation where r ^ x 
above expression simplifies and the 
ward the steady-state (|38|l is Sr{t) 

1 r2[(a-/3)lnr-|-f /3cose]-^{f cos Sin a;}" 
Int [ a+l3+2M. {In x-|-ir(7-ln 2)} _ 

these results tell us that the 2D system evolves loga- 
rithmically slowly toward a space-dependent fluctuating 
steady state. As in the presence of only one zealot, we 
can see that in 2D the magnetization does not exhibit a 
scaling expression between r and t when ^ t ^ 1 0. 
This is due to the logarithmic terms, specific to the two- 
dimensional situation, appearing in H35I) and l|36|) . Nat- 
ural questions arise regarding the spatial distribution of 
"opinions" : What is the spatial voting distribution in the 
steady-state? Which region is characterized by a major- 
ity of positive/negative opinion ? How does the strength 
of a and (3 affect the final spatial opinion distribution ? 

To answer these questions, we use Ea. (|57|l and notice 
that in the limit r ^ 1 and |r — a;| ^ 1, the spatial 



region where Sr{oo) — obeys the equation: 



r — X 



with X = sxp 



n[f3- 



aj3 



(40) 



When a ^ (3, i.e. for x 7^ 1; such an equation can 
be recast into the following form: + ^i^i cos 9 — 

— 0, i.e. the polar equation of a circle C{c,R) 

centered at c — , ^ o X and of radius R = oi — 

x/2sinh(|Q;-i - f3~'^\). This resuh, together with |jSZ|) 
and 1)41) |l . implies that in 2D, for a f3, the agents lo- 
cated on the circle C{c,R) are "neutral" they have zero 
final magnetization as illustrated in Fig. 5. From H37|) 
and H4l)|l we can also conclude that: 

• If X > 1; f3 > a, the agents that are within 
(outside) the disk Int C(c,i?) have a positive (negative) 
magnetization. 

• If X < 1; (3 < a, the agents that are within 
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FIG. 6: The stationary site magnetization S',-(oo) on a (128 x 128) lattice in the presence of six zealots. Three of the zealots 
favor the +1 state and other three the —1 opinion. The picture on the left shows a 3D plot of ^^.(oo) (along the vertical axis) 
and the picture on the right is the corresponding contour plot. The strengths of the positive zealots are 2.0, 1.2,0.8 and the 
strengths of the negative ones are 1.6, 1.4, 1.0. 



(outside) the disk Int C(c,i?) have a negative (positive) 
magnetization. This case is sketched in Figure 2. 

These results show that the majority of the voters, ex- 
cept the ones enclosed in the disk, tend to follow the 
opinion favored by the strongest zealot. The details 
of the neutral region C{c,R) depend nontrivially on all 
the parameters a,/3 and x and, interestingly, the radius 
grows with the difference of the strength of the zealots as 
R (X 1/ sinhu, where u = f3~^ — . Also, R increases 
linearly with the separating distance x. 

• The case a = (5 (including a — (3 — oo), i.e. x = 1, 
is special. In this situation, it follows from (|40|l that the 
region with zero-final magnetization is no longer a closed 
curve but an infinite line given by the equation r = ^^^^^ 
which separates the two-dimensional space into two semi- 
infinite half-planes. 

For the number of zealots n > 2 the analytical calcula- 
tions become very tedious and we illustrate the results of 
a Monte Carlo simulation of the case with six zealots on 
Fig- El The simulation is carried on a (128 x 128) lattice 
and due to the 1/ ln(t) approach to the steady state enor- 
mous sampling times are required. Again when simulat- 
ing the system one has to be careful with the sites on the 
boundaries: if the site lies on the edges then it has only 
three nearest neighbors; and if it is at the corners, then it 
has only two nearest neighbors. The stochastic rules have 
to be slightly modified to account for the boundary sites. 
The geometry of the zealots can be seen from contour 
plot on Fig. El where three of the zealots are positively 
biased and three are biased negatively. The left picture 
on Fig. El shows the average magnetization Sr {oo) on the 
different sites of the lattice. For these particular values 
of the bias of the zealots and their position on the lat- 
tice, in the stationary state, we observe one large region 
of positive on average opinion (a curved central "stripe" 
in Fig. EI) and two smaller disconnected regions of a neg- 
ative opinion (near the left boundary and top right edge 



of Fig. El- 

Regarding the coarsening of the 2D system, we again 
distinguish three situations: (i) when n < 2, there is 
usual coarsening and an infinite domain eventually spans 
the entire system; (ii) when there is a finite number of 
competing zealots large domains still develop but their 
size is limited by the zealots; (iii) when the density of the 
competing zealots is finite in the thermodynamic limit, 
there is no longer coarsening as the formation of large do- 
mains is prevented by the interaction with the numerous 
inhomogeneities. 

To conclude this section, as in ID, we notice that a{\ — 
5*0(00)) = l3{l+Sx[oo)) which implies, with Ea. l(TH|) . that 
the global magnetization evolves, following the strongest 
zealot (a ^ /?), as M(t) ^ i/lnt. As a consequence, the 
time T necessary for the strongest zealot to dominate 
(on average) the whole 2D system is T ^ In L (where 
L — !■ 00). In the symmetric case (a = /?), as explained 
above, the 2D space is exactly separated in two semi- 
infinite half-planes with opposite total magnetization. 



C. Results in 3D 

Above two dimensions, the integrals in Eq. are well 
defined for all values of s and in particular when s — s- 0. 
Therefore, in contrast to what happens in \D and 2D, 
to determine the long-time behavior of the magnetization 
we cannot simply focus on the q ^ Q expansion of 10} . 
This also means that in dimensions d > 3 in the presence 
of n zealots the static magnetization readily follows from 
from Ea. (fTE|l : 

n n 

Sr{^) = q4._,(0)[AA-1(0,M)],-,. (41) 

3 = 1 1=1 
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The three-dimensional lattice Green function 7^(0) has 
been computed very recently by Glasser and Boersma 
[2^ . Using the triplet (0^,6^,0^) of rational numbers 
depending on r, given in Table 2 of Reference and 

the quantity go = (^) (^) C^) = 0.505462... 

(r(z) is Euler's Gamma function), it has been established 
that: 



With (|20|l and (|42|l the exact expression of the three- 
dimensional magnetization in the presence of two zealots 
is explicitly given by: 



t(0) ^ i 



Or .go + Cr -\ 5 



(42) 



5*7.(00) 



aei/.(0) + /3e2t-x(0) + aP {t(0) (ei/o(0) - e2/x(0)) + ir-^{0) (e2/o(0) - eit (0)) } 



-/3)/o(0) + a/?(/2(0) 
I 



(43) 



From now on, for the sake of concreteness, we focus 
on the case where we have two competing zealots, ei = 
—62 = 1, and thus the expression (|43|l becomes Sr{oo) = 

aJ.(0)-ffJ.-^(0)+a/3(/.(0)-t-^(0))(/o(0)+j^(0)) wP arP 

l+{a+f3)io{0)+aP{ii{0)-im) . • WC die 

mainly interested in the large r limit, one can observe 
that /r(0) is just the static solution of the Poisson equa- 
tion Aj.lr{0) = —Sr,Oy which solution in the continuum 
limit is 7^(0) ~ i{r) = (r > 0). This result, ob- 
tained from an "electrostatic reformulation" , is valid on 
the discrete lattice for r ^ 1 [^H. With the help of 
(I43|l . this result allows to compute the 3D stationary lo- 
cal magnetization for r ^ 1 and |r — ccj 1 : 



r 



Co 



where Ci 



2a 



and C2 



2/3 



(44) 



Again, the resem- 



2+ago 2+/3go ' 

blance with electrostatics is striking: the static magne- 
tization is formally the electrostatic potential generated 
by the "charges" Ci at site and C2 at x. As already 
noticed, the difficulty resides in the fact that the charges 
Ci and C2 are a priori unknown and have been computed 
in a self- consistent way (assuming a large enough sepa- 
rating distance a:) , w ith the help of the exact and discrete 
results 631 111. Even though the result (gH) is for- 
mally valid for r ^ a; ^ 1, as explained above, it gives 
already accurate predictions when r ^ 1 and x is finite 
but large enough {e.g. already when x > 6). It is sugges- 
tive that in the limit where a = /3 = 00, the "charges" 
C2 = —Ci = 2/ go. In this case the magnetization in Eq. 
(I44|l can be viewed as the potential of the electric dipole 
of charges ±2/f/o- To make the connection with electro- 
statics even more transparent, it is worthwhile to notice 
that the expression (|44|) can be rewritten using a multi- 
pole expansion. Also, when /3 = 0, we recover 5*^(00) cx 
l/r, as reported in Ref. 0|. In fact, one has \r — x\ 



mials. Thus the expression 1)44(1 can be recast into 



Sr{oo) 



1 

47rr 



00 

Cl+C2 2_,[-) Pmicosd) 



m—O 



(45) 



At this point it is important to mention a major differ- 
ence with the case where only a single zealot is present. 
In the latter situation, as showed in Ref.0, just by taking 
the continuum limit of the equation for the magnetiza- 
tion, one could anticipate that Sr {00) oa r~^ (i.e it has 
only radial dependence) in three dimensions, which is the 
main desired information. In the two-zealot there 
is a competition between the effects of the "charges" Ci 
and C2, we really need to determine Sr {00) through Eqs. 
(|42|l . (|43|) . to obtain the nontrivial spatial dependence of 
the stationary magnetization through H44|l , (|45|l . 

Regarding the dynamical approach toward the steady 
state, it is difficult to study the small s behavior of /r(s) 
and to rigorously obtain the long-time approach toward 
the stationary magnetization. However, it follows from 
Eq. ^ that: 



Sr{t) — S'r(oo) 



1 



2tt (47ri)2 



,(46) 



This result is expected to be accurate in the regime where 
t —^ oo,r 1 and |r — a;| 3> 1. As previously mentioned, 
we can can discuss about the regions with positive or neg- 
ative stationary magnetization. To determine the "neu- 
tral" region (where S'r(oo) = 0) it follows from (|44|l that, 
in the limit where r ^ 1 and |r — a;| ^ 1, one has to 
solve 



with S : 



m + ago) 

a(2 + Pgo) 



(47) 



(r2 



c2 - 2r.a;)-i/2 



cos 9 = ^ and the Pr, 



r ^-^r 

i{cos9 



PjnicosO), where When a ^ P, i.e. for (5^1, the equation can again be 



are the Legendre polyno- recast into the following form: r 



2rx 

'7^ 



cos ( 



= 0. 
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Such an expression is the polar equation of a sphere 
E(C,7?.) centered at C = j^g^x with a radius TZ ~ 
jiZ^Tj- From Ea. H44|) . we can also infer the following: 

• If (5 > 1, i.e. (3 > a, the agents that are within 
(outside) the sphere I](c, i?) have a positive (negative) 
magnetization. 

• If (5 < 1, i.e. (3 < a, the agents that are within 
(outside) the sphere I](C,7?.) have a negative (positive) 
magnetization. 

These results show that majority of the voters, except 
the ones enclosed in the sphere T,{C,TZ), tends to follow 
the opinion favored by the strongest zealot. The details 
of the neutral region S(C, TZ) depend nontrivially on all 
the parameters a,f3 and x. In particular, we notice that 
TZ increases linearly with the separating distance x. 

• The case where a = f3, i.e. 6 = 1, is special because 
thus the "effective charges" are such that \Ci\ = €2- In 
particular, this is the case when a — f3 = 00. It thus 
follows from (|47|) that the region with zero-final magne- 
tization is no longer a surface but an infinite plane, given 
by r = jjf^, that separates the 31? space into two re- 
gions. 

In 3D, a{l ~ 5o(oo)) 7^ /3(1 + S^{oo)) when a ^ /3, 
and thus the global magnetization of the above inho- 
mogeneous voter model evolves linearly with the time: 
M{t) ~ t. This implies that the time T necessary for the 
strongest zealot to dominate (on average) the whole 3D 
system scales as T '--^ L^, where L ^ 00. On the other 
hand, when a — (3, the space is divided in two symmetric 
regions with opposite total magnetization. 

Finally, in the case where both zealots favor the same 
opinion e = ±1, i.e. ei = £2 = e, one has just to modify 
the expressions of "charges" in Eqs 14411 . (|45|l and (|46|l . 
In fact, these results are still valid with Ci — — „^^" and 

' ^ 2+aga 

r - 2£/3 

"-2 - 2+j3go ■ 



V. MONOMER-MONOMER CATALYTIC 
REACTION ON AN INHOMOGENEOUS 
SUBSTRATE 

The other model that we specifically consider in this 
work is the monomer-monomer catalytic reaction. Such a 
process is of considerable interest in many fields of science 
and the technology. In the catalysis the rate of a chem- 
ical reaction is enhanced by the presence of a suitable 
catalytic material, such as the platinum used to catalyze 
the oxidation of carbon monoxide {2CO -f O2 ^ 2CO2) 
\l9j ,23] ■ Because of the numerous and practical implica- 
tions of the catalytic reaction, it is of prime interest to be 
able to model its quantitative and qualitative behavior. 
In general, these processes are described within mean- 
field like approaches where it is assumed that molecules 
are randomly distributed on the substrate 0, |2^ . Spa- 
tial fluctuations and excluded volume constraints are 
thus ignored, despite of the fact that these effects are 
shown to play often a crucial role [2^ . 



In the modeling of catalysis [23, the monomer- 
monomer surface reaction model plays an important part 
at least from a theoretical point of view because the sim- 
plicity of the model allows to address several issues ana- 
lytically, such as the role of the fluctuations % , the 
interfacial roughening [2^, and the diffusion of the ad- 
sorbents ^26]. 

The monomer-monomer catalytic process on an homo- 
geneous substrate is by now well understood and it com- 
prises the following reactions : 

A + — ^ As 

kA 

B + $ — > Bs 

kB 

As + Bs AS t +20. 

The A and B particles impinge upon a substrate with 
rates kA and ks, respectively, adsorb onto vacant sites 
(0) and form a monolayers of adsorbed particles. As and 
Bs. Nearest-neighbor pairs of different adsorbed parti- 
cles, AsBs, react and desorb with rate kr, leaving two 
vacant sites (20) on the substrate. The dynamics on a 
spatially homogeneous substrate is most interesting in 
dimensions d < 2, when kA = ks (otherwise the species 
with the bigger rate will rapidly saturate the substrate). 
In this case there is coarsening on the substrate induced 
by fluctuations and islands of As and Bs particles grow. 
As in Refs. 5, 6], we will consider the reaction- controlled 
limit , where kr <ti kA = ks. This limit turns out to be 
useful from a technical point of view and, most impor- 
tantly, provides qualitatively the same kind of behavior 
as the general case 0, 0, 01 • In the reaction-controlled 
limit, the substrate quickly becomes fully occupied and 
stays covered with As's and Bg's for ever (vacancies are 
immediately refilled). The kinetics of monomer- monomer 
substrate reaction model is therefore a two-state system 
that can be mapped onto the voter model supplemented 
by an infinite-temperature Kawasaki exchange process 

1^. In fact, in the monomer-monomer catalytic re- 
action under consideration. As and ^5 desorb and the 
resulting empty sites are instantaneously refilled either 
by AsBs (no reaction). As As, BsBs (voter dynamics), 
or by BsAs (Kawasaki exchange dynamics at infinite 
temperature). 

Clearly, more realistic situations should include the 
presence of inhomogeneities which could deeply affect the 
properties of the system. In fact, real substrates (in ID 
and 2D) display generally some degrees of spatial het- 
erogeneity which are attributed to imperfections, such 
as dislocations and defects j22| that modify locally the 
interactions on the substrate. In some previous works 
translationally-invariant disordered models for catalysis 
have been considered within mean-field like approaches, 
i.e. rate equations and pair approximation |28l |. In these 
works, it was shown that quenched substrate imperfec- 
tions dramatically affect the dynamics resulting in a re- 
active steady-state. One should emphasize that both 
the physical systems (in this work, the inhomogeneities 
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are not randomly distributed but fixed) and the analytic 
methods (we obtain exact results in arbitrary dimensions, 
while the authors of employed mean-field-like ap- 
proaches) considered here differ from, and are thus com- 
plementary to, those of Rcf. 28] . Also, very recently, an 
equilibrium model for monomer-monomer catalysis on a 
disordered substrate was solved 29]. 

Hereafter we study the static and dynamical effects 
of local inhomogeneities in the monomer-monomer cat- 
alytic reaction-controlled process and show how to take 
advantage of the results obtained for the inhomogeneous 
voter model to infer some exact properties. In fact, we 
consider the genuine nonequilihrium situation where the 
substrate is spatially inhomogeneous, because of the pres- 
ence of a collection of n inhomogeneities located at sites 
{a^}, j = 1, . . . ,n favoring the local adsorption of A's 
or i3's. We show that the inhomogeneities induce spa- 
tially dependent reactive steady-state when n > 1. As 
a substrate, as described in Section 11, we consider an 
hypercubic lattice with (2L -I- 1)** sites and introduce a 
set of parameters e'j taking the values or 1 and con- 
sider, in addition to the usual homogeneous catalytic re- 
action described above, that some inhomogeneities lo- 
cally favor the presence of A via desorption of i?'s (and 
vice versa) through the additional reactions Bs — > As, 

where e' — 1, and As > Bs, where e', — 0. We there- 

fore consider the following homogeneous processes (voter 
+ infinite-temperature Kawasaki dynamics), all occur- 
ring with the same rates 1/2, and local (inhomogeneous) 
reactions at sites and , occurring respectively 

with rates a,- and 



y '■ 



AsBs 
AsBs 
As 



1/2 
7/2^ 



AsAs; AsBs > BsBs; 

1/2 

BsAs; 



Bs\ Bs 



a' , ; e' — 1 



A 



s- 



Here, the bimolecular reactions correspond to the 
voter dynamics supplemented by Kawasaki infinite- 
temperature exchange process, whereas monomolecular 
processes correspond to reactions induced by local inho- 
mogeneities favoring the adsorption of one species. Fol- 
lowing the same steps as in Refs |^ , for this spatially 
inhomogeneous monomer-monomer catalytic process, in 
the thermodynamic limit we obtain the following equa- 
tion of motion for the concentration c^. {t) of ^5 at site r 
of the substrate: 



di 



Cr{t) — ArCr{t) 



J2 



(48) 



Of course, the concentration of Bs at site r is simply 
given by 1 — Cr{t). The resemblance of Equation (|48|l with 
© is striking (the only difference is that here = 0, 1) 
and one can immediately infer the solution of H48|) from 
(| 151 and H2U|) '). In particular, in the thermodynamic limit. 



starting from a system initially completely occupied by 
Bs particles, the Laplace transform of the concentration 
of As reads: 

C'-(^) = lT.^e^'^^'r{s)W-HsAa})hM (49) 

and we get for the time-dependent concentration (ini- 
tially c^(0) = 0): 

-1 />c-\'ioo J 

= ^ / (50) 
In this language, the quantity 

M'{t) = J2 cfe(i) = E - (^)] (51) 

fe 1=1 ^0 

provides the average total number of the As particles on 
the substrate at time t. 

Next, we restrain ourself to physical situations and 
consider in detail the monomer- monomer catalytic reac- 
tion in the presence of one and two inhomogeneities in 
one and two dimensions. 



A. Inhomogeneous monomer-monomer catalytic 
reaction in the presence of one single "defect" 

Here, we consider the case where there is a single in- 
homogeneity at site = 0, with strength ai = a and 
ei = 1. In this case, we simply have Af'^ = 

^ t- J l+aIo(s} 

Therefore, starting from a system initially full of ,65 par- 
ticle {i.e. Cr{0) = 0) we obtain: 



. . ^ 1 air 

Cr[s) = . 

s 1 + aIois) 



(52) 



On the right-hand side of this equation, one recognizes 
immediately the same expression as the Laplace trans- 
form of the magnetization obtained in Ref. Q. From 
previous results, we can immediately infer the long-time 
behavior of the concentration of As particles. 



1. Results in ID 

Following the same steps as in Ref. 0, on a one- 
dimensional substrate we find from H52I) that the long- 
time behavior of the concentration of ^5 reads: 



Crit) 



2/a 



TTt 



(53) 



This result is valid for any < r < cxi. 

When both r — > 00 and t 00, we obtain the following 
simple scaling expression |^: 



Crit) 



(54) 
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2. Results in 2D 

In two dimensions, following results from Ref.Q we ob- 
tain a non-scaling expression for the concentration, with 
very slow time-dependence: 



co{t) - co(oo) ^ - [ — hp— 



47r\ 1 
a I In t ' 



(55) 



where co(oo) = 1. For the other sites, we find that the 
long-time behavior in the regime t 3> ^ 1 is given by 



In r 
InT 



(56) 



As in the one-dimensional case, the stationary concen- 
tration of ^5 corresponds again to a substrate fully cov- 
ered with ^5 particles, i.e. Cripo) = 1. Therefore, the 
presence of a single inhomogeneity favoring locally the 
adsorption of As is enough to completely cover the sub- 
strate with As in spite of the fact that initially only Bs 
particles were present. From the expressions (|5S|) . (j55|) 
and H51|) . we can also compute the total number of As 
particles on the substrate at time i 3> 1. In so doing, 
one obtains M'{t) ^ \/t in the one-dimensional case and 
M'{t) ~ t/lnt in 2D. 



B. Inhomogeneous monomer- monomer catalytic 
reaction in the presence of two defects 



Here, we consider the case where two "competing" in- 
homogeneities are present: one is at site = 0, with 
strength ai = a and ei = 1 and the other at site — x, 
with strength a2 = (3 and 62 = 0. 

In this case, using Eas. (|49|l and Ullfl . we obtain the fol- 
lowing expression for the Laplace transform of the con- 
centration of As at site r, starting from Cr{0) = 0: 



airjs) + af} {ir{s)io{s) - ir-a:is)ia;{s)) 

s[l + ia + P)io{s) + al3{il{s) - il{s))] ' 



(57) 



1. Results in ID 

In one dimension, without loss of generality, we assume 
that the inhomogeneity at site = a;, x = |a;|, is on the 
right side of the origin. 

Proceeding as in section IV. A, we study the static 
and long-time behavior of the concentration of As with 
Crifi) = 0, and distinguish various situations: 

• For sites between the two inhomogeneities, i.e. < 



r < X we get: 

a[l + /3(x-r)] 



Cr(t) 



a 



a + (3 + afix \a (3 + af3x 
[1 + (3{x - r)]{2 - aPx^ /2) 
a + l3 + al3x 



(58) 



• At the right of the origin, when a; < r < 00, we 
obtain: 



Cr{t) 



a 



a-^ 13 + a/3x 



1 - 



1 



2 - af3x^/2 
nt { a + f3 + a(3x 



(59) 



C-r{t) 



• At the left of the origin, when < r < 00, we find: 

a 



a-\- (3 + af3x 



l + /3.-^(.+ (i±Mtit^ 
VTTt L a + fi + afix 



(60) 



• When both t — > cx) and r — > 00, we have: 



Cr{t) 
C-r{t) 



6rfc f I 

a + l3 + al3x \2^/i J 

"(^ + ^") erfcf^^ 
a + l3 + al3x \2y/t J 



(61) 
(62) 



These results show that in the interval between the of 
inhomogeneities, the static concentration profiles varies 
linearly from the origin with a slope —a[3/{a-\-l3-\-a(3x). 
Outside from this interval, the static concentration is uni- 
form on the right and left side of the origin: on the 
right, c,-(oo) = -^^j^^^^, whereas on the left 0^(00) = 
a+i3+a(3x ' ^^'^^ ^ static profile can again be interpreted 
as the solution of a discrete ID electrostatic Poisson 
equation with peculiar and suitable boundary conditions. 
Again, the static concentration is reached according to a 
power-law {cr{t) ~ t^^/^) and with amplitudes depend- 
ing nontrivially on all parameters of the system. At very 
large distances, and long time, the concentration displays 
a scaling form which amplitude depends on which inho- 
mogeneity is the closest. Of course, it is easy to check 
that in the limit a — *■ , as the system is initially full of 
Bs, then Cr{t) = 0. Also, when /3 — 0, we recover the 
expressions and ltM|) . From Eqs (|5T|l and we 
obtain the average number of adsorbed particles which 
evolves (at long-time) as M'(t) ~ ^/t. 

Again, in one dimension we can obtain the sta- 
tionary concentration of adsorbed ^5 particle in the 
completely disordered case, i.e when n is arbitrary 
large just by replacing respectively S^j (00), S'r(oo), by 
(00), Cr(oo), Ej- in the expressions l|77|l - l|55|) . As iUus- 
trated in Fig.01 in this case the stationary concentration 
profile is piecewise. Also, when the number of competing 
inhomogeneities is finite the system coarsens as described 
in Section IV. A. 
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2. Results in 2D 



In two dimensions and at large distance from both 
inhomogeneities, i.e for r 3> 1 and |r — a;| 3> 1, we 
find a non-scaling expression for both static and time- 
dependent concentration of the As particles: 



a/3 
2tt 



In 



2;|r— ic| 



- 7r(7 - In 2) 



a + 13 +^{\nx + Tr{j- In 2)) 



(63) 



and, when x is large enough, co(oo) 

Cx(oo) 



a{l+^ Inx) 
a+l3+^ Inx 



and 



■ In X 



We can notice that in 2D the stationary concentration 
of the ^5 particles is a fluctuating reactive state exhibit- 
ing nontrivial radial and polar dependence. Regarding 
the approach toward the steady state, proceeding as in 
the section IV. B, we obtain: 



• {t) — Cr{ao) 



B'{r,x) 



(64) 



1 ii 1-4. J m -^lna;{ln 7 — x +7 — In 2}+2a In r 

where the amplitude B = — ' „tir, ; — \ — 

exhibits a nontrivial spatial dependence. Again, the re- 
sult (|54|) shows that the stationary concentration profile 
(|63|l is reached logarithmically slowly. Using Eq. H51|) 
we can also notice that the average number of particles 
As adsorbed on the substrate evolves (at long-time) as 
M'{t) ~ i/lni. 

There is a practical interest in understanding the spa- 
tial distribution of adsorbed particles in the steady state 
[3QJ and one can thus ask: What is the region of the 2D 
substrate where one can find more As particles ? 

To answer this question, from Eq. H63|) , we proceed as 
in Section IV. B and, according to Ea. (|40|l . we see that 
when a > /3 (/3 > a), the region richer in As particles 
is outside (within) the disk Int C(c,i?) [defined in Sec- 
tion IV. B], where the concentration of As is Cr(oo) > ^ 
{Cr {00) < i). When a = /3, the 2D substrate is separated 
into two half-planes with concentration of As > 1/2 in 
the region including the origin. 



VI. SUMMARY AND CONCLUSION 

In this work we have shown how to compute some exact 
properties of a class of many-body stochastic systems in 
the presence of an arbitrary number of inhomogeneities 
n, and have specifically focused on the voter model and 
monomer-monomer catalytic reaction (in the reaction- 
controlled limit). We have studied the effects of local 
perturbations of the dynamical rules on the static and 
time-dependent properties of these models by obtaining 
both general (yet formal) and many explicit results in 
the presence of one and two inhomogeneities. In fact, 
the latter situation already displays and covers most of 
the generic features of the models. Namely, when there 



is only one inhomogeneity present, it is responsible for a 
uniform and "unanimous" steady state in low dimensions 
4], while in the presence of competing inhomogeneities 
{n > 1) the steady state is fluctuating and reactive. For 
the sake of concreteness we have mainly focused on the 
amenable case with two inhomogeneities and have shown 
quantitatively how the local interactions deeply affect the 
properties of these systems. Neither the stationary nor 
the time-dependent expression of the order parameters 
are translationally-invariant but exhibit nontrivial radial 
and polar dependence (when d > 1). 

From a sociophysical perspective, in the voter model 
language, this means that a system which tolerates the 
presence of "competing zealots", i.e. which accepts the 
competition between opposite points of view, will never 
reach a unanimous state but always end into a final con- 
figuration where both opinions coexist and fluctuate. Of 
course, such a conclusion seems to be consistent with the 
results of electoral competitions in modern democracies. 

In the presence of competing inhomogeneities {n > 1) 
in low dimensions, subtle coarsening phenomena take 
place in ID and 2D. In fact, the local and competing 
perturbations of the dynamics lead us to distinguish the 
case where the number of inhomogeneities is finite and 
the case where their number is comparable to the size of 
the system. In the former case the system coarsens and 
large domains develop, but their size are typically lim- 
ited by the number of competing inhomogeneities, while 
in the latter case coarsening is prevented by the interac- 
tion with all the numerous inhomogeneities. 

More specifically, in this work we have obtained exact, 
yet formal, expressions of the static and time-dependent 
order parameters (see (|T^ and (j^ l. The main tech- 
nical problem to carry out detailed calculation resides 
in the inversion of the n x n matrix J\f. The case with 
one single inhomogeneity in the voter model was already 
considered in j3| and here we show that such results can 
be translated in the language of the catalysis reaction. 
In particular we have shown that on ID and 2D sub- 
strates, the presence of a single spatial inhomogeneity 
favoring the adsorption of one species, say As, with re- 
spect to the other is sufficient to ensure that eventually 
the substrate will be completely filled with As particles. 
When we have two competing inhomogeneities, favoring 
locally opposite states or the adsorption of particles of 
different species, we have obtained rich behavior. In ID, 
between the two inhomogeneities, the stationary profiles 
of the order parameters vary linearly with the distance 
from the origin (|22|l . (|58|) and then reaches two plateaus 
(I^SJ , 1231 and ifS^ , (EOl . These static profiles are always 
reached algebraically in ID: Sr{t) — S', (oo) ~ At^^^'^ and 
Cr{t) — Cr(oo) ~ A't'^^^^, whcrc the amplitudes A and 
A' depend nontrivially on all parameters of the problem 
and in particular on the separating distance between the 
inhomogeneities [see Eqs (HSJ, (EH) and ((23, (EOH- Far 
away from the inhomogeneities, the order parameters dis- 
play scaling expression of the variable r/^/i [see H25(l and 
(|61|l . (|62|l ]. In one dimension, we have also been able to 
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compute the expression of the stationary magnetization 
in the completely disordered situation where the number 
of zealots is arbitrary large [see Eqs. l(2IJ-l|SSl)]- In two 
dimensions, for n = 2, in the presence of two compet- 
ing inhomogeneities, we have obtained non-uniform and 
nontrivial stationary profiles for the order parameters, in 
agreement with an electrostatic-like reformulation, the 
latter display logarithmic spatial dependence (radial and 
polar) [ 137(1 and H63() ]. The approach toward the reactive 
steady state is very slow: Sr{t) — Sr{oo) ~ B/lnt and 
Cr{t) — Cr(oo) ~ B' / \nt, with amplitudes B and B' de- 
pending again nontrivially on all parameters of the prob- 
lem [see (|39|l and H64|) ]. In 2D, for the inhomogeneous 
voter model, we have also studied the spatial regions with 
positive/negative static magnetization and have shown 
that only within a circle, whose center and radius de- 
pend on the strength of the "zealots" and on the distance 
between the latter, the sign of the magnetization is the 
one favorite by the "weakest" zealot. When both zealots 
have the same strength, there is positive/negative mag- 
netization in half-space. In three dimensions, for n = 2 
and in the continuum limit, we have shown that the sta- 
tionary magnetization of the inhomogeneous voter model 
displays a radial and polar dependence that can be recast 
into a multipole expansion ((44|l . corresponding formally 
to the electrostatic potential generated by two "charges" 
that are determined self-consistently using exact results 
from the discrete lattice system. The connection with 
electrostatics is particularly striking in the limit where 
both zealots have an infinite strength, thus the stationary 
magnetization corresponds to the potential of an electric 
dipole. The approach toward the static magnetization 
follows a power-law: Sr{t) — Srioo) ~ Ct~^/'^ (see (^Hl)- 
Also, in 3D we have studied the spatial regions with pos- 
itive/negative magnetization and have shown that out- 
side from a sphere whose center and radius depend on 
the parameters of the system and varies linearly with 
the distance separating the zealots, the sign of the final 
magnetization is the one favored by the strongest zealot. 

The results obtained from Monte Carlo simulations of 
one and two-dimensional lattices show excellent agree- 
ment with the theoretical results obtained for an infinite 
system. In the presence of multiple {n > 2) competing 
inhomogeneities the calculations in two dimensions be- 
come very tedious and we consider this case by numeri- 



cal simulations which confirm the extremely slow dynam- 
ics and the existence of nontrivial spatial dependence of 
the order parameters. We also would like to point out 
one intriguing and interesting fact about the small time 
behavior of the magnetization of the zealots in the one- 
dimensional case. As it can be extracted from Fig. O 
So{t) and S'x(t), for small t, evolve as a power law with 
an exponent numerically smaller than 0.50. The small 
time behavior of the site magnetization of the usual one- 
dimensional voter model (no inhomogeneities) is linear, 
i.e. Sr{t) — Sr{0) oc t for any site r on the lattice. We 
think it would be interesting to investigate further this 
"anomalous" small-t behavior of the magnetizations of 
the zealots in the one and the two-dimensional cases and 
we plan to do it in our future work. Various general- 
izations of this work could also be investigated. For 
instance, it would be worthwhile to consider that the 
inhomogeneities would not be fixed but spatially dis- 
tributed according to some function 7'({a-'}). In this 
case, one should also average on the quenched disor- 
der (on the samples) and one would have to compute: 
Srit) (X E{a^T^i{a'})Sri{a^},t), whcrc Sraa'},t) is 
the quantity studied in this work for a given set of inho- 
mogeneities at sites {a-'}- In the same manner, it would 
be quite interesting to consider the disordered case where 
the strength of the inhomogeneities would follow a dis- 
tribution function such as V{{aj}) oc Oj^i e~'^"^~"'> 1"^" . 
In this case, one would be interested in the quan- 
tity: Sr{t) = /njdajP({aj})S',.({aJ},{a^},t), where 
5r ({a-'}, {a-'}, t) is the magnetization computed in this 
work for a given set of inhomogeneities at sites {a-'}, with 
strength {a-'}. 
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